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We investigate the Mott transitions in two-orbital Hubbard systems. Applying the dy- 
namical mean field theory and the self-energy functional approach, we discuss the stability 
of itinerant quasi-particle states in each band. It is shown that separate Mott transitions 
occur at different Coulomb interaction strengths in general. On the other hand, if some 
special conditions are satisfied for the interactions, spin and orbital fluctuations are equally 
enhanced at low temperatures, resulting in a single Mott transition. The phase diagrams 
are obtained at zero and finite temperatures. We also address the effect of the hybridization 
between two orbitals, which induces the Kondo-like heavy fermion states in the intermediate 
orbital-selective Mott phase. 

§1. Introduction 

Strongly correlated electron systems with some orbitals have been investigated 
extensively. 1 )' 2 ) In particular, substantial progress in the theoretical understanding 
of the Mott transition in multiorbital systems has been made by dynamical mean- 
field theory (DMFT) 3 )~ 7 ) calculations. 8 )~ 40 ) Among them, the orbital-selective Mott 
transition (OSMT) 41 ) has been one of the most active topics in this context. A 
typical material is the single-layer isovalent ruthenate alloy Ca2- ;c Sr. c Ru04. 42 )- 44 ) 
The end-member Sr2Ru04 is a well-known unconventional superconductor, 45 )' 46 ) 
while Ca2Ru04 is a Mott-insulating 5 = 1 antiferromagnet. 47 )" 49 ) The relevant Ad- 
orbitals belong to the ^-subshell. The planar structure prevents the hybridization 
between orbitals which have even (d xy ) and odd parity (d yz ,d zx ) under the reflection 
z — > —z. The complex evolution between these different end-members has stimulated 
theoretical investigations, 50 )~ 53 ) and among others to the proposal of the OSMT: 
some of the (i-orbitals fall in localized states, while the others provide itinerant 
electrons. The OSMT scenario could explain the experimental observation of a 
localized spin S = 1/2 in the metallic system at x ~ 0.5 in Ca2_ a; Sr :r Ru04, which is 
difficult to obtain from the entirely itinerant description. 41 )' 52 )' 54 ) Another example 
of the OSMT is the compound La n +iNi n 03 n +i. 55 ) It is reported that the OSMT 
occurs in the e 9 -subshell at the critical temperature T c ~ 550K, below which the 
conduction of electrons in the 3d x 2_ y 2 orbital is disturbed by the Hund coupling 
with the localized electrons in the 3d 3z 2_ r 2 orbital. 56 ) These experimental findings 
have stimulated theoretical investigations of the Mott transitions in the multiorbital 
systems. 8 )- 20 )' 41 )' 52 )' 54 ) 

In this paper, we give a brief review of our recent studies 14 )' 15 )' 20 )' 33 )' 36 )' 40 ) on 
the Mott transitions in the two-orbital Hubbard model by means of DMFT and the 
self-energy functional approach (SFA). 57 ) In particular, we focus on the role of the 
orbital degrees of freedom to discuss the stability of the metallic state at zero and 
finite temperatures. The paper is organized as follows. In §2, we introduce the model 



2 



A. Koga, K. Inaba, and N. Kawakami 



Hamiltonian for the two-orbital systems and briefly explain the framework of DMFT 
and SFA. We first treat the Hubbard model with same bandwidths, and elucidate 
how enhanced spin and orbital fluctuations affect the metal-insulator transition in 
§3. In §4, we then consider the system with different bandwidths to discuss in which 
conditions the OSMT occurs. We obtain the phase diagrams at zero and finite 
temperatures. Finally in §5, the effect of the hybridization between the two orbitals 
is investigated to clarify the instability of the intermediate OSM phase. A brief 
summary is given in the last section. 

§2. Model and Methods 

2.1. Two-orbital Hubbard model 

We study the two-orbital Hubbard Hamiltonian, 

H = H + H' (2-1) 
#o = E i Q)c L^ CT + ^ E [ c lia^ + 4 2a c tl a} (2-2) 

<i,j> ia 
aa 

H' = J2 H 'i (2-3) 

i 

H'j = Uy^, ri ia ^n iai + E ~ J '^aa') nii a n i2 a' 

a aa' 

~ 3 E C L C ^ C L C *2a - J' 4lT C !l| C i2tCi2| + cj^cj^C^C^ (2-4) 



where c\ aa {ci acr ) creates (annihilates) an electron with spin a{=],[) and orbital 
a(= 1,2) at the ith site and 

T^iaa — cj acr Qacr* ^ij^ ^ the orbital-dependent nearest- 
neighbor hopping, V the hybridization between two orbitals and U (U') represents 
the intraorbital (interorbital) Coulomb interaction. J, J z and J' represent the spin- 
flip, Ising, and pair-hopping components of the Hund coupling, respectively. We note 
that when the system has the Ising type of anisotropy in the Hund coupling, J = 
J' = 0, the system at low temperatures should exhibit quite different properties from 
the isotropic case. 8 ) -13 )' 27 )' 38 ) In this paper, we deal with the isotropic case and set 
the parameters as J = J z = J', which satisfy the symmetry requirement in the multi- 
orbital systems. It is instructive to note that this generalized model allows us to study 
a wide variety of different models in the same framework. For 1^ = 0, the system is 
reduced to the multi-orbital Hubbard model with the same (t^ = Uj) or distinct 

orbitals. 8 )~ 20 ) On the other hand, for = 0, the system is reduced to a correlated 
electron system coupled to localized electrons, such as the periodic Anderson model 
(J = 0), 58 )~ 62 ) the Kondo lattice model (V = and J < 0) 63 )' 64 ) for heavy-fermion 
systems, and the double exchange model (V = and J > 0) for some transition metal 
oxides. 65 )~ 68 ) For general choices of the parameters, various characteristic properties 
may show up, which continuously bridge these limiting cases. 
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2.2. Method 

2.2.1. dynamical mean- field theory 

To investigate the Hubbard model (2-4), we make use of DMFT, 3 )- 7 ) which 
has successfully been applied to various electron systems such as the single-orbital 
Hubbard model, 69 )~ 82 ) the multiorbital Hubbard model, 8 )~ 40 ) the periodic Ander- 
son model, 83 ) -91 ) the Kondo lattice model, 92 ) -94 ) etc. In the framework of DMFT, 
the lattice model is mapped to an effective impurity model, where local electron 
correlations are taken into account precisely. The Green function for the original 
lattice system is then obtained via self-consistent equations imposed on the impurity 
problem. 

In DMFT for the two-orbital model, the Green function in the lattice system is 
given as, 



G(k,z)- 1 = G (k,zy 1 -V(z), (2-5) 

where (i is the chemical potential, and e a (k) is the bare dispersion relation for the 
a-th orbital. In terms of the density of states (DOS) p(x) rescaled by the bandwidth 
D a , the local Green function is expressed as, 



with 



and 



Gn(z)= I dx P{X) , 



l(*) = / - 
n ( \ f \v p{x)v{z) 

GMz) = J dx ^( Z ,x)aM-v(zr 

G 22 (z) = [ dx P{X) , (2-8) 



where 



£i(z,x) = z + [i- Un - Dix, 
£ 2 (z, x) = Z + fl - E 22 ~ D 2 x, 

v(z) = V + U 12 (z) . (2-9) 

In the following, we use the semicircular DOS, p(x) = ^Vl — x 2 , which corresponds 
to the infinite-coordination Bethe lattice. 

There are various numerical methods to solve the effective impurity problem. 
We note that self-consistent perturbation theories such as the non-crossing approx- 
imation and the iterative perturbation method are not efficient enough to discuss 
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orbital fluctuations in the vicinity of the critical point. In this paper, we use numer- 
ical techniques such as the exact diagonalization (ED) 71 ) and the quantum Monte 
Carlo (QMC) simulations 95 ) as an impurity solver at zero and finite temperatures. In 
this connection, we note that the Hund coupling is a crucial parameter that should 
control the nature of the Mott transition in the multiorbital systems. 9 ) It is thus 
important to carefully analyze the effect of the Hund coupling in the framework of 
QMC. For this purpose, we use the QMC algorithm proposed by Sakai et al., 37 ) 
where the Hund coupling is represented in terms of discrete auxiliary fields. When 
we solve the effective impurity model by means of QMC method, we use the Trotter 
time slices At = (TL) -1 < 1/6, where T is the temperature and L is the Trotter 
number. 

2.2.2. self-energy functional approach 

We also make use of a similar but slightly different method, SFA, 57 ) to determine 
the phase diagram at finite temperatures. This SFA, which is based on the Luttinger- 
Ward variational method, 96 ) allows us to deal with finite-temperature properties 
of the multi-orbital system efficiently, 20 )' 40 ) where standard DMFT with numerical 
methods may encounter some difficulties in practical calculations when the number 
of orbitals increases. 

In SFA, we utilize the fact that the Luttinger-Ward functional does not depend 
on the detail of the Hamiltonian TLq as far as the interaction term 7i' is unchanged. 57 ) 
This enables us to introduce a proper reference system having the same interaction 
term. One of the simplest reference systems is given by the following Hamiltonian, 

^ref = Yli ^ref ' 



n 



(0 
ref 



CWT 



e ( V c 

e 0a L icwr c *' 



+ e 



W o (0t o « 



,0 A 



,0 

0.(7 



l«t, 



+ ^,(2-10) 



where o iaa {a-atr) creates (annihilates) an electron with spin a and orbital a, which is 
connected to the ith site in the original lattice. This approximation may be regarded 
as a finite-temperature extension of the two-site DMFT. 76 ) In the following, we fix 
the parameters eo a = and e a = fi to investigate the zero and finite temperature 
properties at half filling. We determine the parameters v a variationally so as to 
minimize the grand potential, dfi/dv a = (a = 1, 2), which gives a proper reference 
system within the given form of the Hamiltonian (2-10). 

§3. Mott transition in the degenerate two-orbital model 

We begin with the two-orbital Hubbard model with same bandwidths (tff = Uj) 
at half filling. 21 ) -40 ) In this and next sections, we neglect the hybridization term by 
putting V = 0. We first discuss the Fermi-liquid properties in the metallic phase 
when the Coulomb interactions U and U' are varied independently in the absence of 
the Hund coupling (J = 0). 33 )' 36 ) The effect of the Hund coupling will be mentioned 
in the end of this section. 40 ) 
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3.1. zero-temperature properties 

To investigate the Mott transition at zero temperature, we make use of the ED 
method as an impurity solver, where the fictitious temperature 71 ) (3 allows us to 
solve the self-consistent equation G\ oc = Gi mp numerically. In this paper, we use 
the fictitious temperature (5{> 50) and the number of sites for the impurity model 
is set as N = 6 to converge the DMFT iterations. We note that a careful scaling 
analysis for the fictitious temperature and the number of sites is needed only when 
the system is near the metal-insulator transition points. To discuss the stability of 
the Fermi liquid state, we define the quasi-particle weight for ath band as 

-l 



1 - lmU a (u! n ) 



Mr. 



(3-1) 



n=0 



where uj n [= (2n + 1)tt/(3] is the Matsubara frequency. 

In Fig. 1, the quasi-particle weight calculated is shown as a function of the 
interorbital interaction U' for several different values of intraorbital interaction U. 
For U = 0, Z decreases monotonically with increasing U', and a metal-insulator 




Fig. 1. The quasi-particle weight Z as a function of the interorbital Coulomb interaction U' . The 
data are obtained by the ED (N = 6) within DMFT. 

transition occurs around U' c ~ 0.9. On the other hand, when U ^ 0, there ap- 
pears nonmonotonic behavior in Z; it once increases with the increase of U', has the 
maximum value in the vicinity of V ~ U, and finally vanishes at the Mott transi- 
tion point. It is somehow unexpected that the maximum structure appears around 
U ~ U', which is more enhanced for larger U and U' . We will see below that this is 
related to enhanced orbital fluctuations. By repeating similar calculations system- 
atically, we end up with the phase diagram at zero temperature (Fig. 2), where the 
contour plot of the quasi-particle weight is shown explicitly. At U' = 0, where the 
system is reduced to the single-orbital Hubbard model, we find that the intraorbital 
Coulomb interaction U triggers the metal-insulator transition at U c = 2.9 ~ 3.0. 
This critical value obtained by the ED of a small system (N = 6) is in good agree- 
ment with other numerical calculations such as the numerical renormalization group 
(U c = 2.94), 73 )' 81 ) the ED (U c = 2.93), 77 ) the linearized DMFT (U c = 3), 75 )' 76 ) the 
dynamical density-matrix renormalization group (U c = 3.07). 80 )' 82 ) 
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Fig. 2. The zero-temperature phase diagram for J = 0. The contour plot of the quasi-particle 
weight Z is shown: the bold-dashed line represents the phase boundary of the metal-insulator 
transition, which is obtained by estimating the values of U and U' that give Z — 0. The solid 
circles are the transition points obtained by the linearized DMFT. 75 ' 

There are some notable features in the phase diagram. First, the value of Z 
is not so sensitive to U' for a given U (> U'), except for the region U ~ U'. In 
particular, the phase boundary indicated by the dashed line in the upper side of the 
figure is almost flat for the small U region. The second point is that when U ~ U' the 
metallic phase is stabilized up to fairly large Coulomb interactions, and it becomes 
unstable, once the parameters are away from the condition U = U'. The latter 
tendency is more conspicuous in the regime of strong correlations. 24 )' 33 ) 

3.2. finite-temperature properties 

To observe such characteristic properties around U = U' in more detail, we com- 
pute the physical quantities at finite temperatures by exploiting the QMC method 
as an impurity solver. 95 ) In Fig. 3, we show the DOS deduced by applying the 
maximum entropy method (MEM) 97 )~") to the Monte Carlo data. 36 -* In the case 



U'=0.0 U'=3.0 U'=4.0 




CO 

Fig. 3. The DOS for the two-orbital Hubbard model (U = 3.0). The data are for the inverse 
temperature T _1 = 1, 2, 4, 8 and 16 from the top to the bottom. 

(U,U') = (3.0,0.0), the system is in the insulating phase close to the transition 
point (see Fig.2). Nevertheless we can clearly observe the formation of the Hubbard 
gap with the decrease of temperature T. In the presence of U', the sharp quasi- 
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particle peak is developed around the Fermi level at low temperatures (see the case 
of (U, U') = (3.0, 3.0)). This implies that orbital fluctuations induced by U' drive the 
system to the metallic phase. Further increase in the interorbital interaction sup- 
presses spin fluctuations, leading the system to another type of the Mott insulator 
in the region of U' > U. 

To characterize the nature of the spin and orbital fluctuations around the Mott 
transition, we investigate the temperature dependence of the local spin and orbital 
susceptibilities, which are defined as, 

Xs = / drdn^-n^jin^-n^T)}) 
Jo 

Xo= / dT<{m(0)-n 2 (0)}{ni(T)-n 2 (T)}), (3-2) 
J o 

where (5 = T _1 , n a = n QiCT , n a = n a>a and r is imaginary time. We show the 
results obtained by QMC simulations within DMFT in Fig. 4. Let us first look at 




Fig. 4. The local spin and orbital susceptibilities as a function of the temperature. 



Figs. 4 (a) and (b) for U' = (equivalent to the single-orbital model). Since the in- 
troduction of the intraorbital interaction U makes the quasi-particle peak narrower, 
it increases the spin susceptibility Xs at low temperatures. On the other hand, the 
formation of the Hubbard gap suppresses not only the charge susceptibility but also 
the orbital susceptibility Xo- As seen from Figs. 4 (c) and (d), quite different behav- 
ior appears in the susceptibilities, when the interorbital interaction U' is increased. 
Namely, the spin susceptibility is suppressed, while the orbital susceptibility is en- 
hanced at low temperatures. This tendency holds for larger U' beyond the condition 
U' ~ U. Therefore in the metallic phase close to the Mott insulator in the region 
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of U > U'(U < U'), spin (orbital) fluctuations are enhanced whereas orbital (spin) 
fluctuations are suppressed with the decrease of the temperature. These analyses 
clarify why the metallic phase is particularly stable along the line U = U'. Around 
this line, spin and orbital fluctuations are almost equally enhanced, and this subtle 
balance is efficient to stabilize the metallic phase. When the ratio of interactions 
deviates from this condition, the system prefers either of the two Mott insulating 
phases. 

3.3. phase diagram at finite temperatures 

QMC simulations are not powerful enough to determine the phase diagram at 
finite temperatures. To overcome this difficulty, we make use of a complementary 
method, SEA, 40 ) which allows us to discuss the Mott transition at finite temperatures. 
The obtained phase diagram is shown in Fig. 5, where not only the case of J = 
but also J = 0.1U are studied under the condition U = U' + 2J. 



0.1 

T 



0.05 



4 6 8 10 

Fig. 5. The finite temperature phase diagram for the two-orbital Hubbard model for J = and 
J = 0.1(7 under the condition U = U' + 2 J. 

It should be noted that the Mott transition at finite temperatures is of first order, 
similarly to the single-orbital Hubbard case. 3 -* The critical value for the transition 
point, U c , is determined by comparing the grand potential for each phase. We find 
the region in the phase diagram where the metallic and Mott insulating phases 
coexist. Starting from the metallic phase at low temperatures, the increase of the 
Coulomb interaction triggers the Mott transition to the insulating phase at U C 2, 
where we can observe the discontinuity in the quasi-particle weight Z. On the other 
hand, the Mott transition occurs at U c i when the interaction decreases. The phase 
boundaries U c , U c \ and U& merge to the critical temperature T c , where the second 
order transition occurs. 

Note that upon introducing the Hund coupling J, the phase boundaries are 
shifted to the weak-interaction region, and therefore the metallic state gets unstable 
for large U. Also, the coexistence region surrounded by the first order transitions 
shrinks as J increases. This tendency reflects the fact that the metallic state is stabi- 
lized by enhanced orbital fluctuations around U = U': the Hund coupling suppresses 
such orbital fluctuations, and stabilizes the Mott-insulating phase. Another remark- 
able point is that the Mott transition becomes of first-order even at zero temperature 
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in the presence of J^ 24 )' 35 )- 38 ) since the subtle balance realized at T = in the case 
of U = U' is not kept anymore for finite J. There is another claim that the second 
order transition could be possible for certain choices of the parameters at T = 0, 38 ) 
so that more detailed discussions may be necessary to draw a definite conclusion on 
this problem. 

§4. Orbital-selective Mott transitions 

In the previous section, we investigated the two-orbital Hubbard model with 
same bandwidths, and found that the metallic state is stabilized up to fairly large 
Coulomb interactions around U = U' , which is shown to be caused by the en- 
hanced spin and orbital fluctuations. We now wish to see what will happen if we 
consider the system with different bandwidths, which may be important in real ma- 
terials such as Ca 2 _ x Sr x Ru04 42 ) and La n+1 Ni n 03 n+ i. 55 )' 56 ) In the following, we will 
demonstrate that the enhanced spin and orbital fluctuations again play a key role 
in controlling the nature of the Mott transitions even for the system with different 
bandwidths. 14 )' 15 )' 20 ) 

4.1. zero-temperature properties 

Let us start with the quasi-particle weight calculated by DMFT with the ED 
method 14 ) and see the stability of the metallic phase at zero temperature. Here, 
we include the Hund coupling explicitly under the constraint U = U' + 2J. The 
quasi-particle weights obtained with fixed ratios U' /U and J/U are shown in Fig. 6 
for half-filled bands. We first study the case of U = U' and J = with bandwidths 




Fig. 6. The quasi-particle weights Z\ and Z2 at half filling as a function of U for D\ — 1.0 and 
D 2 = 2.0: (a) U' /U = 1.0 (J = 0) and (b) U' /U = 0.5 (J/U = 0.25). Open (closed) circles 
represent the results for orbital a = 1(2) obtained by combining DMFT the ED as an impurity 
solver for the N = 6 small cluster. Solid triangles represent the Mott-transition points obtained 
by the two-site DMFT method. Insets show the same plots for D\ = 1.0 and D2 = 5.0. 

D\ = 1.0 and D2 = 2.0 [Fig. 6 (a)]. When the Coulomb interaction is switched 
on, the quasi-particle weights Z\ and Z2 decrease from unity in slightly different 
manners reflecting the different bandwidths. A strong reduction in the quasi-particle 
weight appears initially for the narrower band. However, as the system approaches 
the Mott transition, the quasi-particle weights merge again, exhibiting very similar 
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[/-dependence, and eventually vanish at the same critical value. This behavior is 
explained as follows. For small interactions, the quasi-particle weight depends on the 
effective Coulomb interactions U/D a which are different for two bands of different 
width D a , and give distinct behavior of Z\ and Zi- However, in the vicinity of the 
Mott transition, the effect of the bare bandwidth is diminished due to the strong 
renormalization of the effective quasi-particle bandwidth, allowing Z\ and Zi to 
vanish together. 100 )' 101 ) 

The introduction of a finite Hund coupling J makes U 7^ U', and causes qualita- 
tively different behavior [Fig. 6 (b)]. As U increases with the fixed ratio U'/U = 0.5, 
the quasi-particle weights decrease differently and vanish at different critical points: 
U c i ~ 2.6 for Z\ and U C 2 ~ 3.5 for Z2- We thus have an intermediate phase with one 
orbital localized and the other itinerant, which may be referred to as the intermedi- 
ate phase. The analogous behavior is observed for other choices of the bandwidths, 
if J takes a finite value [see the inset of Fig. 6 (b)]. These results certainly suggest 
the existence of the OSMT with U C 2 > U c ±. 

We have repeated similar DMFT calculations for various choices of the param- 
eters to determine the ground-state phase diagram, which is shown in Fig. 7. The 

8 
6 
U4 
2 

°0 2 4 6 8 

IT 

Fig. 7. The zero-temperature phase diagram for two-orbital Hubbard model for Di = 1 and D2 = 2. 
In the phase (I) [phase (II)], both bands are in metallic (insulating) state. The phase (III) is 
induced by the OSMT, where the metallic state coexists with the Mott insulating state. Since 
we consider the ferromagnetic Hund coupling, J > 0, the relevant region in the diagram is 
U > U'. 

phase diagram has some remarkable features. First, the metallic phase (I) is stabi- 
lized up to fairly large Coulomb interaction U when U — > V (small J) . Here the Mott 
transitions merge to a single transition. As mentioned above, this behavior reflects 
the high symmetry in the case of U = U' (J = 0) with six degenerate two-electron 
onsite configurations: four spin configurations with one electron in each orbital and 
two spin singlets with both electrons in one of the two orbitals. Away from the sym- 
metric limit, i.e. U > U' (2 J = U — U') orbital fluctuations are suppressed and the 
spin sector is reduced by the Hund coupling to three onsite spin triplet components 
as the lowest multiplet for two-electron states. In this case, we encounter two types 
of the Mott transitions having different critical points. In between the two tran- 
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sitions we find the intermediate OSM phase (III) with one band localized and the 
other itinerant. Within our DMFT scheme we have confirmed that various choices 
of bandwidths give rise to the qualitatively same structure of the phase diagram as 
shown in Fig. 7 (see also the discussions in Summary). 

4.2. finite-temperature properties 

To clarify how enhanced orbital fluctuations around U = U' affect the nature 
of Mott transitions, let us study the temperature-dependent orbital susceptibility, 

3 



2 
Xo 
1 



°0 1 2 3 4 5 
T 

Fig. 8. The orbital susceptibility as a function of T for V = 0. Open (solid) symbols represent the 
results in the case U = U' and J = (U'/U = 3/4 and J/U = 1/8) and dashed lines those for 
the non-interacting case. 

which is shown in Fig. 8. Here, we have used the new algorithm of the QMC 
simulations proposed by Sakai et a/. 37 ) to correctly take into account the effect of 
Hund coupling including the exchange and the pair hopping terms. 15 ) 

We can see several characteristic properties in Fig. 8. In the non-interacting 
system, the orbital susceptibility increases with decreasing temperature, and reaches 
a constant value at zero temperature. When the interactions are turned on (U'/U = 
3/4 and J/U = 1/8), the orbital susceptibility is suppressed at low temperatures, 
which implies that electrons in each band become independent. Eventually for U > 
U& ~ 3, one of the orbitals is localized, so that orbital fluctuations are suppressed 
completely, giving Xo = at T = 0. On the other hand, quite different behavior 
appears around U' = U. In this case, the orbital susceptibility is increased with 
the increase of interactions even at low temperatures. Comparing this result with 
the phase diagram in Fig. 7, we can say that the enhanced orbital fluctuations are 
relevant for stabilizing the metallic phase in the strong correlation regime. While 
such behavior has been observed for models with two equivalent orbitals (previous 
section), it appears even in the systems with nonequivalent bands. 36 ) 

In Fig. 9, the effective Curie constant % S T is shown as a function of the tempera- 
ture. We first look at the case of U'/U = 3/4 and J/U = 1/8. At high temperatures, 
all the spin configurations are equally populated, so that the effective Curie constant 
has the value 1/2 for each orbital in our units, giving XsT ~ 1. For weak electron 
correlations (U = 1), the system is in the metallic phase, so that the Pauli paramag- 
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netic behavior appears, resulting in XsT — ► as T — ► 0. We can see that the increase 
of the interactions enhances the spin susceptibility at low temperatures, as a result 
of the progressive trend to localize the electrons. The effective Curie constant is 
X S T = 2 when a free spin is realized in each orbital, while % S T = 8/3 when a triplet 
S = 1 state is realized due to the Hund coupling. It is seen that the Curie constant 
increases beyond these values with the increase of the interactions (U = 3,4). This 
means that ferromagnetic correlations due to the Hund coupling appear here. 

For U' = U, we can confirm that not only orbital but also spin fluctuations 
are enhanced in the presence of the interactions, see Fig. 9 (b). Accordingly, both 
spin and orbital susceptibilities increase at low temperatures, forming heavy-fermion 
states as far as the system is in the metallic phase. Note that for U = 6, at which the 
system is close to the Mott transition point, the spin susceptibility is enhanced with 
the effective Curie constant XsT ~ 4/3 down to very low temperatures [Fig. 9 (b)]. 
The value of 4/3 originates from two additional configurations of doubly-occupied 
orbital besides four magnetic configurations, which are all degenerate at the metal- 
insulator transition point. Although not clearly seen in the temperature range shown, 
the Curie constant XsT should vanish at zero temperature for U = U' = 6, since the 
system is still in the metallic phase, as seen from Fig. 7. 

To see how the above spin and orbital characteristics affect one-electron prop- 
erties, we show the DOS for each orbital in Fig. 10, which is computed by the 
MEM. 97 )"") When the interactions increase along the line U'/U = 3/4 and J/U = 
1/8, the OSMT should occur. Such tendency indeed appears at low temperatures in 
Fig. 10(a). Although both orbitals are in metallic states down to low temperatures 
(T = 1/6) for U = 1, the OSMT seems to occur for U = 2; one of the bands develops 
the Mott Hubbard gap, while the other band still remains metallic. At a first glance, 
this result seems slightly different from the ground-state phase diagram shown in 
Fig. 7, where the system is in the phase (I) even at U = 2. However, this deviation 
is naturally understood if we take into account the fact that for U = 2, the narrower 
band is already in a highly correlated metallic state, so that the sharp quasi-particle 
peak immediately disappears as the temperature increases. This explains the be- 
havior observed in the DOS at T = 1/6. For U = 3, both bands are insulating 
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Fig. 10. Solid (dashed) lines represent the DOS for the orbital a = 1 (a = 2) when (Di,D2) = 
(1.0, 2.0). The data are for the temperatures T = 2, 1, 1/2 and 1/6 from the top to the bottom. 

at T = 1/6 (the system is near the boundary between the phases (II) and (III) at 
T = 0). For U' = U, the qualitatively different behavior appears in Fig. 10. In this 
case, quasi-particle peaks are developed in both bands as the interactions increase, 
and the system still remains metallic even at U = U' = 3. As mentioned above, all 
these features, which are contrasted to the case of U' / U, are caused by equally 
enhanced spin and orbital fluctuations around U = U'. 

4.3. phase diagram at finite temperatures 

Having studied the spin and orbital properties, we now obtain the phase diagram 
at finite temperatures in the general case U ^ U' and J / 0. Since each Mott tran- 
sition at zero temperature is similar to that for the single-orbital Hubbard model, 3 ) 
we naturally expect that the transition should be of first order at finite temperatures 
around each critical point. 11 ) In fact, we find the hysteresis in the physical quanti- 
ties. For example, we show the entropy per site in Fig. 11. At T = 0.002, when U 
increases, the metallic state (I) disappears around U C 2 ~ 2.36 where the first-order 
Mott transition occurs to the intermediate phase (III), which is accompanied by the 
jump in the curve of entropy. On the other hand, as U decreases, the intermediate 
phase (III) is stabilized down to U c \ ~ 2.24. The first-order critical point U c ~ 2.33 
for T = 0.002 is estimated by comparing the grand potential for each phase. 20 ) The 
phase diagram thus obtained by SFA is shown in Fig. 12. It is seen that the two 
coexistence regions appear around U ~ 2.4 and U ~ 3.3. The phase boundaries U c i, 
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Fig. 11. The entropy as a function of U at T = 0.002, J = 0.2517. 
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Fig. 12. The finite-temperature phase diagram for J = 0.25C/. There are two coexistence phases, 
which correspond to the triangular-shaped regions: the metallic phase and the intermediate 
phase coexist in the left region and the intermediate phase and insulating phase coexist in the 
right region. 



U C 2 and U c merge at the critical temperature T c for each transition. We note that 
similar phase diagram was obtained by Liebsch by means of DMFT with the ED 
method. 11 ^ Our SFA treatment elucidates further interesting properties such as the 
crossover behavior among the competing phases (I), (II) and (III). To see this point 
more clearly, we show the detailed results for the entropy and the specific heat in Fig. 
13. There exists a double-step structure in the curve of the entropy, which is similar 
to that found at T = 0, where the residual entropy S/L = 0,log2 and log 3 ap- 
pears in the phase (I) , (III) and (II) , respectively. Such anomalies are observed more 
clearly in the specific heat. It is remarkable that the crossover behavior among three 
phases is clearly seen even at high temperatures. Therefore, the intermediate phase 
(III) is well defined even at higher temperatures above the critical temperatures. 

§5. Effect of hybridization between distinct orbitals 

In the present treatment of the model with DMFT, the intermediate phase (III) 
is unstable against certain perturbations. There are several mechanisms that can 
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Fig. 13. (a) Entropy S/L and (b) Specific heat C/L as a function of U in the crossover region, 
J = 0.2517. 



stabilize this phase. A possible mechanism, which may play an important role in 
real materials, is the hybridization between the two distinct orbitals. The effect 
of the hybridization is indeed important, e.g. for the compound Ca2- x Sr x Ru04, 42 ) 
where the hybridization between {a, (3} and 7 orbitals is induced by the tilting of 
Ru06 octahedra in the region of Ca-doping 0.2 < x < O.5. 43 ) An interesting point is 
that the hybridization effect could be closely related to the reported heavy fermion 
behavior. 42 ^' 44 ) The above interesting aspect naturally motivates us to study the 
hybridization effect between the localized and itinerant electrons in the intermediate 
phase (III). Here, we take into account the effect of hybridization in each phase to 
study the instability of the OSMT. 15 )' 18 ) 

We study the general case with U' 7^ U and J / in the presence of the 
hybridization V. In Fig. 14, we show the DOS calculated by QMC and MEM for 
several different values of V. We begin with the case of weak interaction, U = 1, 
where the metallic states are realized in both orbitals at V = 0. Although the 
introduction of small V does not change the ground state properties, further increase 
in V splits the DOS, signaling the formation of the band insulator, where all kinds 
of excitations have the gap. On the other hand, different behavior appears when the 
interactions are increased up to U = 2 and 3. In these cases, the system at V = 
is in the intermediate or Mott-insulating phase at T = 1/6. It is seen that the DOS 
around the Fermi level increases as V increases. At U = 2, the intermediate state 
is first changed to the metallic state, where the quasi-particle peaks emerge in both 
orbitals (V = 0.75, 1.0). For fairly large V, the system falls into the renormalized 
band insulator (V = 1.5). In the case of U = 3, the hybridization first drives the 
Mott-insulating state to the intermediate one, as seen at V = 0.75, which is followed 
by two successive transitions. 

These characteristic properties are also observed in the charge, spin and orbital 
susceptibilities at low temperatures (Fig. 15). For weak interactions (U = 1), the 
charge susceptibility \c monotonically decreases with the increase of V. When elec- 
tron correlations become strong, nonmonotonic behavior appears in Xc'- the charge 
fluctuations, which are suppressed at V = 0, are somewhat recovered by the hy- 
bridization. For large V, Xc is suppressed again since the system becomes a band 
insulator. It is seen that the orbital susceptibility Xo shows nonmonotonic behavior 
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Fig. 14. Solid (dashed) lines represent the DOS for the orbital a = 1 (a = 2) when (Di,Z)2) = 
(1.0, 2.0) at T = 1/6 with the fixed parameters of U' /U = 3/4 and J/U = 1/8. The data are 
plotted for V = 0.0, 0.25, 0.5, 0.75, 1.0, 1.25 and 1.5 from top to bottom. 



similar to the charge susceptibility, the origin of which is essentially the same as in 
Xc- In contrast, the spin susceptibility decreases with the increase of V irrespective 
of the strength of the interactions. As is the case for V = 0, the effective spin is 
enhanced by ferromagnetic fluctuations due to the Hund coupling in the insulating 
and intermediate phases. When the hybridization is introduced in these phases, the 
ferromagnetic fluctuations are suppressed, resulting in the monotonic decrease of the 
effective Curie constant. 

We can thus say that the introduction of appropriate hybridization gives rise 
to heavy-fermion metallic behavior. Such tendency is observed more clearly in an 
extreme choice of the bandwidths, (Di,!^) = (1.0,10.0), as shown in Fig. 16. 
Since the system is in the intermediate phase at V = 0.0, the narrower band shows 
localized-electron properties [Fig. 16 (b)] in the background of the nearly free bands. 
This double structure in the DOS yields two peaks in the temperature-dependent 
Curie constant, as shown in Fig. 16 (a). The localized state plays a role of the /- 
state in the Anderson lattice model, 102 ) so that heavy-fermion quasi-particles appear 
around the Fermi level for finite V, which are essentially the same as those observed 
in Fig. 14. 

Finally, we make some comments on the phase diagram at T = 0. Although we 
have not yet studied low-temperature properties in the presence of V, we can give 
some qualitative arguments on the phase diagram expected at zero temperature. As 
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Fig. 15. The charge, spin and orbital susceptibilities as a lunction of V at T = 1/6. 




Fig. 16. (a) Effective Curie constant as a function of T and (b) DOS in the narrower band (a — 1) 
at T = 1/4 for an extreme choice of the bandwidths, (Di,D 2 ) = (1.0, 10.0). The DOS for the 
wider band is not shown here. The other parameters are U = 4.0, U' = 3.0 and J = 0.5. 



shown above, the metallic phase (I) is not so sensitive to V in the weak-V regime. 
This is also the case for the insulating phase (II), where a triplet state (S = 1) formed 
by the Hund coupling is stable against a weak hybridization. There appears a subtle 
situation in the intermediate phase (III). The intermediate phase exhibits Kondo- 
like heavy fermion behavior at low temperatures in the presence of V. However, 
we are now dealing with the half-filled case, so that this Kondo-like metallic phase 
should acquire a Kondo-insulating gap due to commensurability at zero temperature. 
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Therefore, the intermediate phase (III) should be changed into the Kondo-insulator 
with a tiny excitation gap in the presence of V at T = 0. Accordingly, the sharp 
transition between the phases (II) and (III) at V = still remains in the weakly 
hybridized case. 18 ) This is different from the situation in the periodic Anderson 
model with the Coulomb interactions for conduction electrons as well as localized / 
electrons, where the spin-singlet insulating phase is always realized in its parameter 
range. 87 )' 90 ) Consequently we end up with the schematic phase diagram (Fig. 17) for 
the two-orbital model with the hybridization between the distinct orbitals. Recall 
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Fig. 17. The schematic phase diagram for the model with finite V at T = 0. Solid lines represent 
the phase boundaries between the metallic and insulating phases. Dashed line indicates the 
phase boundary between the Mott insulator and the Kondo insulator. 

that the phase has three regions on the line V = 0: (I) metallic, (II) insulating 
and (III) intermediate phases. The metallic phase (I) for small U is simply driven 
to the band-insulator (IV) beyond a certain critical value of V. The intermediate 
phase (III) at V = is changed to the Kondo-insulator (V) in the presence of any 
finite V. As V increases, this insulating state first undergoes a phase transition 
to the metallic phase (I), and then enters the band-insulator (IV). On the other 
hand, the Mott insulating phase (II) first shows a transition to the Kondo insulator 
(V) , which is further driven to the metallic phase (I) and then to the band-insulating 
phase (IV). Note that at finite temperatures above the Kondo-insulating gap, we can 
observe Kondo-type heavy fermion behavior in the intermediate phase with finite V. 

§6. Summary 

We have investigated the Mott transitions in the two-orbital Hubbard model by 
means of DMFT and SFA. It has been clarified that orbital fluctuations enhanced in 
the special condition U ~ V and J ~ play a key role in stabilizing the correlated 
metallic state. In particular, this characteristic property gives rise to nontrivial 
effects on the Mott transition, which indeed control whether the OSMT is realized 
in the two-orbital model with different bandwidths. 

We have demonstrated the above facts by taking the system with the different 
bandwidths (Di,!^) = (1.0,2.0) as an example. In the special case with U = U' 
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and J = 0, the metallic state is stabilized up to fairly large interactions, resulting 
in a single Mott transition. The resulting single transition is nontrivial, since we 
are concerned with the system with different bandwidths. On the other hand, for 
more general cases with U 7^ V and J / 0, the Hund coupling suppresses orbital 
fluctuations, giving rise to the OSMT. We have confirmed these results by computing 
various quantities at zero and finite temperatures. 

Recently, it was reported that when the ratio of the bandwidths is quite large, 
the OSMT occurs even in the case of U = U' and J = 0. 17),18) This result implies 
that the detailed structure of the phase diagram in some extreme cases depends on 
the parameters chosen, and is not completely categorized in the scheme discussed in 
this paper. This naturally motivates us to establish the detailed phase diagrams at 
zero and finite temperatures by incorporating various effects such as the magnetic 
field, the crystalline electric field, 13 ) the lattice structure, etc. 

In this paper, we have restricted our discussions to the normal metallic and Mott 
insulating phases without any long-range order. It is an important open problem to 
take into account such instabilities to various ordered phases, which is now under 
consideration. 
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